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7^ , Abstract: We study the constraints imposed by neutrino oscillation experiments on the 

minimal extension of the Standard Model that can explain neutrino masses, which requires 
the addition of just two singlet Weyl fermions. The most general renormalizable couplings 
of this model imply generically four massive neutrino mass eigenstates while one remains 
massless: it is therefore a minimal 3 + 2 model. The possibility to account for the confirmed 
solar, atmospheric and long-baseline oscillations, together with the LSND/MiniBooNE and 

10 reactor anomalies is addressed. We find that the minimal model can fit oscillation data 

including the anomalies better than the standard 2>v model and similarly to the 3 + 2 

^—^ phenomenological models, even though the number of free parameters is much smaller 

^ than in the latter. Accounting for the anomalies in the minimal model favours a normal 

^ hierarchy of the light states and requires a large reactor angle, in agreement with recent 

^ measurements. Our analysis of the model employs a new parametrization of seesaw models 

that extends the Casas-Ibarra one to regimes where higher order corrections in the light- 
heavy mixings are significant. 
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1. Introduction 

The simplest extension of the Standard Model that can account for the observed neutrino 
masses involves the addition of two extra singlet Weyl fermions. Such extension encom- 
passes very different possibilities depending on the hierarchy of scales between the Dirac 
and Majorana masses. If the latter are zero the model corresponds to two massive Dirac 
neutrinos; in the opposite limit when they are much larger than the electroweak scale, we 
get the Type-I seesaw model [1, 2, 3, 4]. More exotic possibilities are the so-called direct 
or inverse seesaws [5, 6] when the hierarchy is such that a global lepton number symmetry 
is approximately preserved [7]. Finally, if the Majorana mass scale is low but not zero one 
gets the so-called mini-seesaw models [8, 9], that can be considered the minimal models 
with extra light sterile neutrinos [10]. 

The most general renormalizable Lagrangian including two singlet Weyl fermions, com- 
patible with the SM gauge symmetries is given by: 



where the Yukawa matrix y is a 3 x 2 matrix and Mjj is a symmetric matrix of dimension 
two. We can choose a basis where the neutrino mass matrix takes the form 



where Mjsf = Diag(Mi,M2) is diagonal, while my is a 3 x 2 complex matrix. It can be 
shown that generically the spectrum contains four massive states and one massless state. 
The mixing is described in terms of four independent angles and three CP violating phases. 
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Model 


# Am^ 


7^ Angles 


# Phases 


3 V 


2 


3 


1 


3+2 MM 


4 


4 


3 


3+2 PM 


4 


9 


5 



Table 1: Number of independent mass differences, mixing angles and CP phases (only those that 
can enter oscillations of active states) in the standard iv scenario, the 3+2 phenomcnological model 
(PM) and the 3+2 minimal model (MM). 



It should be stressed that the number of free parameters in this model is much smaller 
than in the models usually referred to as the 3 + 2 model that has been recently revived 
as viable possibility to explain the LSND/MiniBooNE and reactor anomalies [11, 12, 13, 
14, 15]. In those models a generic mass matrix of dimension 5 is considered and therefore 
the ultraviolet completion involves necessarily a more complex extension of the Standard 
Model. For example, one that involves heavier Majorana neutrinos some of which are 
integrated out, or several more Weyl species that pair up to make more than three Dirac 
neutrinos. In the following we will refer to the 3 + 2 models usually considered in the 
literature as 3 + 2 phenomcnological model (PM), while the model with just two Weyl 
singlets is called 3 + 2 minimal model (MM). Table 1 summarizes the number of physical 
parameters that can affect oscillations of active neutrinos in the 3+2 PM and MM models 
compared with the standard ?>v model. 

In Ref. [10] , the constraints of neutrino oscillation data on the 3 + 1 MM and 3 + 2 
MM were studied. The former is completely excluded . For the latter, the parameter space 
in the case of degenerate Majorana masses (Mi = M2 = M) was fully explored, and it 
was concluded that only the regions M < 10^9(10"^°) eV or M > 0.6(1.6) eV for the 
normal (inverted) hierarchy, respectively, are presently allowed by data. Furthermore, the 
degenerate limit with M ~ eV did not improve the fits including the LSND/MiniBooNE 
or new SBL reactor data in contrast with the 3 + 1 PM model. On the other hand, a 
qualitative analysis showed that the non-degenerate case with Mi 7^ M2 ~ 0(eV) could 
actually improve those fits, similarly to the 3 + 2 PM. The main objective of this paper 
is to quantify to what extent this is true. Analyses of the non-degenerate case have also 
been considered in [16, 17]. These studies have used the Casas-Ibarra parametrization that 
assumes an approximate decoupling of light and heavy sectors. For heavy masses in the 
0(eV) range, we find however that the Casas-Ibarra limit is not sufficiently precise and we 
have to use a more general parametrization. 

We will not consider in this paper the limit M2 ^ Mi (3+1+1 minimal model), which 
corresponds to the proposal of [18], recently reconsidered in [16, 19]. 

The structure of the paper is as follows. In Sec. 2 we introduce a new parametrization 
of seesaw-type models that extends the one by Casas-Ibarra to a situation where the de- 
coupling of light and heavy sectors is not large enough. In Sec. 3 we present the results of 
the global fit to the 3+2 MM and compare it with that of the 3+2 PM [13, 14], and with 
the standard model. In Sec. 4 we show the present constraints in the 3+2 MM on the 
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physical combinations that will be measured in future neutrino oscillation experiments. In 
Sec. 5 we conclude. 



2. Parametrization of seesaw models beyond Casas-Ibarra 

There are many possible parametrizations of the mass matrix. A good choice will usually 
be one that satisfies two properties: 1) it contains all independent parameters and no 
more, 2) it is convenient to impose existing constrains. Since two mass splittings have been 
determined with high accuracy, ie. the solar and atmospheric one, it makes sense to use a 
parametrization that uses as parameters those physical quantities. 

In the case when my <^ M^r, a convenient parametrization that does precisely this is 
the one first introduced by Casas-Ibarra [20], which exploits the approximate decoupling of 
the light and heavy sectors, using as parameters not only the masses but also the mixings 
that have already been measured. 

If Mi ~ 0{eV), however, the Casas-Ibarra parametrization (which involves a pertur- 
bative expansion in my/M) is not sufficiently precise. In fact, next order corrections in the 
my/M expansion are significant. In the following we present an alternative parametriza- 
tion, inspired in Casas-Ibarra, that coincides with it in the proper limit, but that does 
not assume any expansion in my or and therefore can be used in the full parameter 
space. For an alternative parametrization with similar properties see [21]. 

There should be a 5 x 5 unitary matrix such that 



We can reduce one dimension by projecting out the zero mode, that is we easily find a 3 x 3 
unitary matrix Uq such that Ai^ of the form 



where now uid has zeros in the first row. A possible choice is to take the first row of Uq 
of the form eijk(my)ii{my)2k properly normalized. The remaining rows can be anything 
as long as it is unitary. 

In this way we reduce one dimension of the problem and find a 4 x 4 unitary matrix, 
V, satisfying 



Mu = U* Diag{0, 1712,7713, Ml, M2) 





(2.2) 




1 
V 



) 



t 



(2.3) 



Defining rrii = Diag{m2,m3) and M/^ = Diag{Mi, M2) and 





in block form, eq. (2.3) implies 



A*7niA^ + B*MhB^ = 0^1 = -M~^^^{B*y\A*)miA^B^)~ 




1/2 



(2.5) 
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Therefore a la Casas- Ibarra we can define a 2 x 2 orthogonal matrix 



-1 

h 



which allows to rewrite B in terms of A: 

The unitarity of V implies: 

AA^ + BB'^ = (i + m]''^RHl^^Rm]''^) = A~^{A^y 



The key point is now to use the polar decomposition of A 

A = WH, = H, W^W = I, 
in the previous equation to get 

i7-2 = / + m]''^ R'^ M^^ Rm]''^ . 
Using unitarity, the blocks C and D can also be rewritten as: 

= -A'^BD^ 



and using the polar decomposition of D = WH we find 



(2.6) 



(2.7) 



(2.^ 



(2.9) 



(2.10) 



(2.11) 



(2.12) 



It can be shown that W is unphysical because it can be reabsorbed in a rotation of the 
sterile states, so we can fix it to the identity. 

The full unitary matrix can then be written as 



U 



Uaa Uas 
Usa Uss 



(2.13) 



with 



UpMNS 



1 
H 



Usa=i[0 HM^ ''^Rm 



1 Uas 



iUpMNS 

H. 



(2.14) 



where Upmns is a generic 3x3 unitary matrix resulting from a combination of Uq and W, 
while H and H are defined in eqs. (2.10) and (2.12). 

More concretely, the physical parameters are chosen to be: the four mass eigenstates, 
three angles and two phases in Upmns- 
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and one complex angle in R: 



f cos(6'45 + Z745) sin(6'45 + Z745) \ ^2 
\^ - sin(6'45 + ^745) cos(6'45 + ^745) J ' 

We note that, even though Upmns is unitary and we have used the same notation as for 
the standard unitary mixing matrix in the model, the mixing in the 3x3 hght sector 
depends also on H, which contains the expected non-unitarity effects. To leading-order 
in the seesaw limit, H H I and in this case the mixing matrix simplifies to the 
Casas-Ibarra parametrization as expected. 

( UpMNS iUpMNsrnJ'^R)M^^^'^\ 

where Upmns is the standard ?>v mixing matrix. 

For the NH case, this is the standard parametrization. However for the IH, since 
the order of the mass eigenstates is unchanged, the parametrization of Upmns is not the 
standard one: it differs from it by the permutation 123 — )• 231, which in particular implies 
that the role of the solar angle is played by the angle ^13 and the role of the reactor angle 
is controlled by 612 (with 612 ~ 7r/2). 

It is straightforward to extend this general parametrization to the hr = 3 case. 

The physical range of the parameters can be chosen to be 9ij G [0, 7r/2], 6,a G [0,27r] 
and 745 G (—00, 00). In contrast with Casas-Ibarra, no divergence occurs for large 745 since 
the mixing matrix converges exponentially to a finite result when this parameter grows. 

The model can obviously be matched to a 3-1-2 phenomenological model, with Am^^ = 
M^,Am|^ = M|, that are free parameters, and with the 5x5 unitary matrix given by 
eqs. (2.13) and (2.14). The elements Uei, Ue5, U^^ and J/^s are well known functions of the 
chosen parameters, and are strongly correlated among themselves and with the masses. In 
Ref. [10, 16, 17], the expression for these quantities in various limits was given using the 
Casas-Ibarra parametrization. 

The expressions for the more general parametrization beyond Casas-Ibarra limit are 
not as simple. We can derive, however, some robust results in various limits. For the NH 
in the limit of vanishing solar mass splitting (i.e. m2 = 0), the heavy-light mixings simplify 
significantly to 



UeA = i 



"^e^(-^)5i3sinz45 + O(vWiW^) , 



Ue5 = n/-^e*^""'^5l3COSZ45 + 0(vWM) 



Uf,4: = "^^^ e'"ci3S23 sinz45 + 0{\/m2/Mi) , 

C^mS = iJ^^^^^^^e'^cviS2i cos Z45 + O(vWM) • (2.18) 
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where 2:45 = ^45 — 1745 and 

X = M1M2 + m3(Mi - M2)/2cos(26'45) + ms^Mi + M2)/2cosh(2745) . (2.19) 

It is easy to check that the quantity cp^^ = aicg{U*^Ue5U^4U*^), which drives the CP vi- 
olation in the fj, — e probabihty in the LSND/MiniBooNE range, is exactly zero in this 
limit. 

We also see that, generically, the electron couplings are suppressed by ^13. These 
couplings are those that get modified more significantly by the effects of the solar splitting. 
In the Casas-Ibarra limit we find the first correction to be 



Uei^i[ ^^Si3SmZ45 + J^Ci3Si2COSZ45 



'e5 



— e ^Sl3COSZ45 - W— Ci3Si2SmZ45 ) . (2.20) 



CP violation in (p4^ starts therefore at 0(^711,2513) or 0(7712) ■ Both terms can be comparable 
in size. Based on these analytical results, it was argued that if the hierarchy is normal (NH), 
the heavy-light mixings would be too small to reach the best fit points found in [13, 14]. 
We will see that this is not the case. The NH works thanks in part to the recent indication 
of a large reactor angle from Double CHOOZ, Daya Bay and RENO, and in part to the 
fact that actually corrections to the Casas-Ibarra limit are significant. 

The ^ couplings, on the other hand, are not small for the NH. In order to keep them 
small enough, there must be also a partial cancellation with higher order corrections in 
eq. (2.18). It is generic that if such a cancellation works for fi, it becomes an enhancement 
in the r channel and viceversa. 

For the IH case, the limit of vanishing solar splitting is more complicated and the 
phase (j)45 does not vanish in general. In fact, even the limit of vanishing solar splitting 
and reactor angle (i.e. m2 = 7713 ~ \/ Am^^^, 612 = 7r/2) is non-trivial in this case. In the 
Casas-Ibarra limit we get 



Uei - '^^Si3SinZ45 + C13 COS 2:45 



TT ■ 1^3 f i(a-S) 

y M~ \ ~ '^^^ ^™ 

Uf,4 ^ ^y^^"^'^23 (^e'(""'')ci3SinZ45 - Si3COS2;45 

Uf,5 ^ ^y§-e''^^23 (e'("-'')ci3 COS Z45 + S13 sin 245) , (2.21) 

where in the IH case S13 corresponds to the solar angle and is therefore not small. The 
mixings are in the right ballpark in this case for Mj ~ eV [10]. 

3. Oscillation data and the 3 + 2 minimal model 

As we have seen, the 3 + 2 minimal model can be represented as a 3 + 2 PM, as the previous 
parametrization shows explicitely. On the other hand, the parameters in the mixing matrix 
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3 + 2 PM 


|Am2,|(eV2) 


|AmiJ(eV2) 


\UeA\ 




\U^.4\ 


\U^.5\ 


045 


KMS[13] 


0.47 


0.87 


0.128 


0.138 


0.165 


0.148 


1.64 vr 


GL[14] 


0.9 


1.61 


0.13 


0.13 


0.14 


0.078 


1.51 TT 



Table 2: Best fit values of tlic licavy mass splittings in the 3 + 2 phenomenological model fits from 
[13, 14]. 

are strongly correlated. The active-sterile mixing, for example, generically depends on all 
the angles and masses, light and heavy. 

The analysis of oscillation data in the context of the 3+2 PM [13] is based on the 
sensible assumption that most of the parameters can be assigned to one of two groups: the 
ones that contribute mostly to long-baseline (LBL) ocillation data (including KamLAND, 
MINOS, CHOOZ, T2K, solar and atmospheric data), such as the two lighter mass split- 
tings, and the ones that contribute mostly to short-baseline (SBL) data (short baseline 
reactors, LSND/MiniBooNE), such as the heavier mass splittings. Only the combinations 
de = |C/e4p + |C^e5p and djj, = |C//x4p + If^A^sP significantly constrained by both [22]. 

The LBL constraints on these combinations are added to the Xsbl puhs and from the 
combined x^i the heavy masses and the active-sterile mixings are extracted. For more 
details see Ref. [22, 13]. 

The decoupling of the LBL and SBL analyses does not hold so cleanly in the 3 + 2 
MM, because the number of independent parameters in the mixing matrix is significantly 
smaller and correlations are large. A global fit would require to deal with all parameters: 
however, this is too time consuming and therefore some approximations are needed. The 
parametrization proposed above reflects the approximate decoupling of the heavy sector 
when Mi ^ rrij that can obviously be exploited. In particular, for NH, the parameters: 
012,(^23, m2,m3, which dominantly drive the solar and atmospheric oscillation, are very 
well constrained by the LBL data^. In order to fix them we will therefore minimize only 
xIbl- The data samples we include as LBL are: KamLAND [23], MINOS-CC [24], T2K 
[25]. In our SBL data sample we follow the same procedure as in [13] and include: LSND 
[26], MiniBooNE [27, 28], KARMEN [29], NOMAD [30], CDHS [31], as well as the reactor 
data Bugey3 [32], Bugey4 [33], ROVNO [34], Krasnoyarsk [35], Gosgen [36], CHOOZ [37] 
and Palo Verde [38] . Oscillations at these shorter baselines do depend also on the previous 
parameters (see eqs. (2.18-2.21)), but SBL data is less constraining, because of larger 
uncertainties both theoretically (dependence on more free parameters) and experimentally. 

On the contrary, the heavy masses Mi and M2 are only relevant at short baselines. 
Therefore they can be determined by minimizing Xsbl- This has been already done in 
the context of the 3 + 2 PM fits [13, 14], and relatively small regions were found for these 
masses. We will therefore use the best fit points of those fits to fix Mi, M2. The results of 
the best fit points we consider are summarized in Table 2. 

The remaining five parameters (6*13(^12)1 ^5 a, 6*45, 745) for NH(III) can affect signifi- 
cantly both LBL, SBL oscillations, and obviously also the reactor oscillation recently dis- 

^For IH the angle relevant for solar oscillations in 613 instead of 612. 
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covered by Double CHOOZ [39], Daya Bay [40] and RENO [41]. We will use a combined 

= xIbL + XSBL + XoayaBay + XreNO ^ t^^em. 

Summarizing, the global is computed in a five-dimensional grid, after minimizing 
Xlbl o'^^r the parameters m2, m^, ^13(^12), ^23 for NH(IH) and fixing Mi, M2 to their best 
fit values in Table 2. The remaining free parameters are allowed to vary within their full 
physical ranges. 

In Fig. 1 we show the results, for the NH and IH, of the one-dimensional projection of 
the total s-s a function of the complex angle 745 , that is the parameter better constrained 
in both cases. The projections on the other parameters are not particularly illuminating, 
because there are strong correlations. In the next section we will consider instead the 
constraints on several directly measurable combinations of the mixing matrix elements. 




-3-2-10123-2 -1 1 2 

745 745 



Figure 1: Minimum in the 3+2 MM as a function of the parameter 745 for NH (left) and IH 
(right), and for A/j = ^/Amj^, M2 = \/Amf^ fixed to the best fit values of the 3+2 PM KMS fit 
(up) and to the GL fit (down). The dashed line corresponds to the minimum of the standard 
3v model. The band is the 2a band of the 3+2 PM KMS/GL best fit. 

As shown in Fig. 1, the 3+2 MM outperforms the standard 31/ scenario, for both NH 
and IH, and for the two choices of Mi and M2 in Table 2. We get a better fit for the KMS 
choice of Mi and M2 for NH, while the GL values of Mi and M2 give a slightly better fit 
in the IH case. The shift in the minimum x^ ^or the two best cases is found to be 

Ax' = xLn(3 + 2 MM) - xLn(3z^) ~ 28 (21) NH(IH) (3.1) 

which is significant given that the difference in the number of free parameters in both 
models is five (two of which are phases). However, the improvement comes solely from the 
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3 + 2 MM 


|Am2j(eV2) 


|AmiJ(eV2) 


\Ue4\ 


\Ue5\ 




\U^.5\ 


4>45 


NH 


0.47 


0.87 


0.149 


0.127 


0.112 


0.127 


1.8 vr 


IH 


0.9 


1.61 


0.139 


0.122 


0.138 


0.107 


1.4 vr 



Table 3: Best fit values of the heavy mass splittings in the 3 + 2 MM for NH and IH. The masses 
are not fitted, are fixed to the best fit points of the phenomenological fits. 

SBL data sample, while the fit to the LBL data deteriorates. This is reflecting the well- 
known tension between appearance and disappearance data. We find that MINOS CC 
sample imposes very significant constraints. The reason is that the MINOS near detector 
is too near and no significant effect of the heavier splittings is visible there, but the far 
detector CC sample is depleted at high energies due to the presence of the heavier states. 
The data do not show any indication of such depletion. 

As an indication, we find that xLn(LBL + SBL) - xLn(LBL) - xLn(SBL) ~ 16(19) 
for NH (IH) for the KMS choice of Mi, Mg. 

Compared to the 3+2 PM, the fit of the minimal model is necessarily worse, since it 
is equivalent to a 3 + 2 PM with constraints. On the other hand Fig. 1 shows that the fit 
does not degrade very significantly, especially in the case of NH and the KMS choice of 
Ml and M2. This is non trivial because the 3+2 MM is significantly more predictive. The 
difference in the number of physical parameters of the 3+2 PM and 3+2 MM is 7 (5 angles 
and 2 phases). 

In Figs. 2 we show the 3a limits on the heavy-light mixing elements \Uai\, i = 4, 5, a = 
e, /i, T in the 3+2 MM. The active-sterile couplings for the electron and muon sectors 
are in the same ballpark as those found in the phenomenological fits. The r mixings, 
which are undetermined in the phenomenological fits, are very constrained in this case and 
turn out to be rather large for the NH. In Fig. 3 we show the limits on the SBL phase 
4>45 = aTg{U*4Ue5Uf^4U*^). The values of the electron and muon mixing with the heavy 
states obtained for the 3+2 MM at the best fit point for NH (KMS masses) and IH (GL 
masses) are shown in Table 3. We note that the corrections to the Casas-Ibarra limit are 
rather significant, up to 30% in some cases. 

One might worry that significant constraints can be obtained from the neutral current 
measurement in MINOS, which we have not included. We have checked that the value of 
the parameter fs of [42], which measures the sterile fraction of the muon disappearance, is 
around 0.2 at the best fit, so it is well within the 90% CL of the MINOS bound. Future 
searches for r appearance oscillations however could be useful to further constrain this 
mode [16]. One would need an experiment able to search for appearance of r's, ie with 
sufficient high energy, and this would require an intermediate baseline of L ~ O(10)km. 
At the oscillation peak the P^r probability would be above 15%. At OPERA one would 
expect a constant contribution in the oscillation probability of ~ 2.5%, which is probably 
too small to be detectable. 
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0.0 0.1 0.2 0.3 0.4 0.0 0.1 0.2 0.3 0.4 

|t/„4l \Uai\ 



Figure 2: Left: Scr ranges on the plane (|C/q4|,|?7q5|) for a — e (solid), /i (dashed) and r (dotted) 
for the NH, with Mi and M2 fixed to the KMS values of Table 2. Right: same for the IH and with 
the GL values of Mi and M2. The symbols: circle (e) and square {^) correspond to the best fit 
points of the 3+2 PM fits from Table 2. 




4. Phenomenological implications for LBL reactors and accelerators 

The LBL reactor and T2K/N0VA measurements aim at the subleading atmospheric oscil- 
lation in the electron channel. In this section we review what to expect in the 3+2 MM 
model, since the previous fits constrain most of these observables. In a more general PM 
model, there is more freedom and some of these predictions can be relaxed. 

In reactor and accelerator experiments the far detector is located at a baseline compa- 
rable to solar or atmospheric oscillations and therefore, being the heavier mass splittings 
significantly larger than the atmospheric one, are in the fast-oscillation regime. In appendix 
A, we give the general oscillation probabilities in this regime. Often these experiments also 



- 10 - 



have near detectors to monitor the neutrino flux and further control other systematic er- 
rors. While the near dectector in reactor experiments are typically also rather far, so that 
the heavier splittings are still in the fast oscillation regime, this is not the case with ac- 
celerator experiments, where the near detectors are at a distance where either oscillations 
have not taken place or they are about to. 

The electron disappearance probability measured by reactors, in the approximation of 
neglecting the solar mass splitting, can be written as 



PI AT 

eel reactor ~ ^^e. 



1 - sm"' ^^^^ 



(4.1) 



with 



Nee = 1-2 {\Ue4\^ + \Ue5\^ - \Ue4\^ - \Ue5\^ " | C/e4| ' | f/e5 1') , 

A..Nee.iTr^!^:^:^\:^rO: (4.2) 



\4\Uel\\\Ue3\' + \Ue2\^). IH 



(4.3) 



Reactor experiments such as Daya Bay and RENO use near detectors to fix the normaliza- 
tion, therefore they are sensitive to the quantity Age- Note that the near detectors in both 
experiments are located sufficiently far to observe averaged oscillations for the two large 
mass splittings . 

The projection of the function of Aee is shown in the top plots of Fig. 4 for 

the NH and IH. It is interesting to see that the 3a limits are more restrictive than the 
corresponding 3a limits of the recent Daya Bay measurement, shown by the shaded region. 

The so-called golden-channel oscillation probability that will be measured in the T2K 
far detector (FD) and other future experiments is given, in the same approximation of one 
scale dominance (atmospheric), by 



(±4> + ^"^""^ ) + 0{Asolar), 



(4.4) 



where 



N^e = 2 {\UeAnUf,i\^ + IC/eslVMsl + K ( f/e4 C/;4 C/^g C//,5 ) ) 

1 4 Ueiu;,{u:,u,3 + u:,u,2). IH 



(4.5) 
(4.6) 



The lb refers to neutrinos and antineutrinos respectively. As explained in the appendix 
there is CP violation even in the one scale dominance limit, e.g A21 — 0. The effect 
is however suppressed by the ammount of violation of unitarity of the light-sector (see 
appendix A) , as expected on general grounds [43] . 
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Figure 4: Minimum as a function of the physical combinations A^e, 2|A 
NH(left) and IH (right). The heavy masses are fixed to the KMS values of Table 2 for NH and to 
the GL values for IH. The band in the top plots corresponds to the recent Daya Bay Sa limit. 



The minimum ^ a function 2|^^e| (which in the standard 3iy scenario is roughly 
Age) and A'^e (which would be zero in the case) are shown in the second and third 



- 12 - 



row of Fig. 4. The constraint on the phase (j) is shown in Fig. 5. As expected, these 
non-conventional CP violating effects are small (we can see that ~ vr for NH and IH). 

The T2K experiment also has a near detector. However the baseline of 280 m is 
actually near the oscillation length corresponding to the mass splittings of Table 2, so 
the approximation of averaged oscillations is not good in this case. A careful analysis 
would require to fit simultaneously the near and far detector oscillation probabilities. The 
appearance probability in the near detector can be approximated by 



P ~ 



4|C/e4|'|C/M4l'sin2 



AE 



+ 4|C/e5|'|C//.5|'sin2 



S\UeAU;^U:^U^,^\ sin ( I sm 



AE 



AE 



cos 



Af|L 
AE 
(M|- 



Ml)L 



AE 



=F <P45 



(4.7) 



while the disappearance one, that can also be measured, is 

|2 \TT |2\ I \TT |2 „• 2 ^^iL 



P^^ ~ l-4(l-|[/^4 

- 4|?7^4|VA.5|'sin^ 



f/M5|') (|i7^4|'sin2 ^ + |[/^5l'sin2 -^2 



AE 



2,^^ ,2 ^^^2 { M^-Ml)L 
AE 



(4.8) 



This is also the good approximation for LSND and MiniBooNE. The corresponding mixing 
elements and the phase have been presented in the previous section. 
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Figure 5: a function of the non-conventional phase for NH (left) and IH (right). The 

heavy masses are fixed to the KMS choice of Table 2 for NH and to the GL choice for IH. 



In Fig. 7 we show the P^e and probabilities for the best fit of the 3+2 MM (NH) 
at a baseline of L = 280m and in the range of energies of the T2K beam. Clearly the T2K 
near detector could provide very valuable information to further constrain the 3+2 MM. 

Finally, in the longer term, the goal in future LBL experiments is to achieve sufficient 
precision to measure also the CP violation driven by the solar and atmospheric interference, 
ie. the only one present in the standard 3z^ scenario. In this case, we cannot consider the 
limit of vanishing solar mass splitting. A simultaneous expansion in Am^^^^^ and the small 
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Figure 7: Appearance probability /ie (left) and fi disappearance (right) at a baseline L = 280m 
(T2K near dectector) for the 3+2 MM best fit values for NH (sohd line) and the 3+2 PM best fit 
values of [13] (dashed line). 



mixings in eq.(A.l) gives the following approximate result 

r> \ AT \ A \ ■ ( ^atmL\ . / , , , ^atm.L 
^/.e|T2K-FD - " I^M^I Sm ( ^ I Sm ( ±(/. + ^ 

, In I (IS- , ^atrriLx . ^solarL . Aatm.L 

+ -Due COS ±Ocff H — sm sm — - — 



COS iCOcff 

+ 0{Al,l,„eAsolar), (4.9) 

with 

\B \,±^s^. ^ h U^^^UU^sU,;, NH 

' \ 8 [/ei[/*i[/*3C/^3e*^. IH ^ - ' 

(4.11) 

The minimum as a function of \Befi\ is also shown in the last row of Fig. 4. The 
constraint on the effective phase dgg (that plays the same role as the standard CP phase in 
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the Sv model) is shown in Figs. 6. We see that the range of values of 5es in the 3+2 MM 
is poorly constrained: whilst a non-vanishing d^s is expected for NH, all values are allowed 
for IH. 

5. Conclusions and Outlook 

We have studied a 3+2 minimal model with light sterile neutrinos, which corresponds 
to enlarging the Standard Model with two extra singlet Weyl fermions (often referred 
to as mini-seesaw models), and is arguably the simplest extension that can explain the 
observed neutrino masses. The neutrino spectrum contains generically four massive states 
and therefore can be represented as a 3+2 phenomenological model, but with constraints. 
Mixings and masses are correlated, as in seesaw models. 

As it is well known, models with extra light sterile neutrinos in the 0(eV) range do 
improve the global fits to neutrino oscillation data, because they can accommodate the 
LSND signal (recent analyses can be found in [13, 14, 15]) and the new reactor anomaly, 
but there remains a significant tension in the /i sector between positive appearance signals 
and negative disappearance ones. In the minimal model, it is possible to fix the same 
physical spectrum as found in the 3+2 phenomenological fits but, once the spectrum is 
fixed, the mixing parameters are severely constrained: only 4 angles and 3 phases can be 
further tuned, while in the phenomenological 3+2 model, 9 angles and 5 phases enter in 
oscillation probabilities. The minimal model is therefore much more predictive. 

We have evaluated to what extent the minimal model can fit global oscillation data, 
in spite of these constraints. We have presented a useful parametrization of seesaw mod- 
els inspired by the Casas- Ibarra parametrization [20], but one that does not assume any 
expansion in the ratio of light and heavy masses. It reduces to Casas-Ibarra in the appro- 
priate limit, but can be safely used even if corrections of higher order in the ratio of light 
to heavy masses are significant, as happens in our case. We believe this parametrization 
can be useful in other contexts. 

The results of the global fits show that the 3+2 minimal model outperforms the 3z^ 
model in global fits by about a shift of 28 (20) units of for NH (IH). Previous estimates 
based in Casas-Ibarra parametrization [10, 16, 17] indicated that the electron mixings to 
the heavy states would be too small for the NH. We arrive to a different conclusion, due in 
part to the large value of the recently measured reactor angle and to the fact that higher 
order corrections in the heavy- light mixings are significant. The minimal model for the NH 
predicts however significantly larger tau mixings to the heavy states. 

Compared to the 3+2 phenomenological fits, we find that the minimal models, in spite 
of the constraints, give similarly good fits, especially in the case of the NH. 

Experimental evidence of extra light sterile neutrinos is slowing accumulating from 
oscillations experiments and cosmology. Unfortunately, the significant signal from LSND 
has not yet been confirmed at a similar level of confidence by any other experiment. If 
such a signal would be confirmed in the future, it might be explained by just adding two 
Weyl singlet fermions to the Standard Model. 
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A. Appendix 



We gather here some useful expressions for the oscillation probabilities in 3+2 models. 

For values oiE/ L such that the heavier splittings are in the averaged-oscillation regime, 
the oscillation probability is well approximated by the following expression 



Pap = Sapil - 2eaa) + 2 (|ea/3p - \Ua4U^4U*^Ui35\ COS (f) 45) 



+ 8\Ua2U*^2Ua3Ul33\ COS ±023 



A21L . AaiL 



sm ■ 



sm 



+ MUasU^s^-aisl sin ( ±03 + — ^ ] sm 



Casl sin I ±02 ± 



A21L 



sm 



2 

A21L 



(A.l) 



with 

3 

= Sai3 - ^ U^i^^i <Pij = arg{f7aiC/^jC/*jC/^j}, 0i = arg{C/aif7^je„/3}. (A. 2) 
1=1 



The e terms measure the violation of unitarity in the light 123 sector and are therefore the 
terms that are non-standard. The probabilities are recovered in the limit e, Ua4, Uab — ^ 0. 

As expected there is CP violation even in the one scale dominance limit, e.g A21 — 
0. The effect is however suppressed in light-sector non- unitarity, as expected on general 
grounds [43]. The corresponding CP violating terms (last two terms in eq. (A.l)) are indeed 
suppressed in 0(e). 

Let us rewrite Eq. (A.l) for a ^ [3 expanding over the solar oscillation frequency. For 
the NH case we obtain: 
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p. 



yNH _ 

aP — 




(A.3) 



where 



and the phases (pi and have been defined above. As expected, the appearance probabihty 
in the IH case can be obtained just doing 123 —t- 231 in the above equations: 



Of course, the same symmetry between NH and IH is shown by the disappearance (and 
appearance) probabiUties presented in Sec. 4. In this way we are expressing the probabihties 
as a function of the solar and atmospheric mass differences since the transformation 123 —)• 
231 corresponds to the fohowing mapping: 



Finahy, note that no expansion over e^^ has been considered here. 
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